Bound states and persistent currents in topological insulator rings 
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We analyze theoretically the bound state spectrum of an Aharonov Bohm (AB) ring in a two- 
dimensional topological insulator using the four-band model of HgTe-quantum wells as a concrete 
example. We calculate analytically the circular helical edge states and their spectrum as well as 
the bound states evolving out of the bulk spectrum as a function of the applied magnetic flux and 
dimension of the ring. We also analyze the spin-dependent persistent currents, which can be used to 
measure the spin of single electrons. We further take into account the Rashba spin-orbit interaction 
which mixes the spin states and derive its effect on the ring spectrum. The flux tunability of the ring 
states allows for coherent mixing of the edge- and the spin degrees of freedom of bound electrons 
which could be exploited for quantum information processing in topological insulator rings. 
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I. INTRODUCTION 

Mesoscopic rings offer an ideal testbed for the investi- 
gation of geometrical phase effects in quantum mechanics 
like the well known Aharonov-Bohm (AB) effect 1 . Ex- 
perimental observation of the energy spectra of quantum 
rings, showing clear signatures of AB periodicity, have 
been obtained by means of optical spectroscopy 2 * 3 . The 
predicted persistent currents^ have been experimentally 
measured, also on a single ring, in various experiments 
in metals^— and semiconductors-. Phase effects, such 
as the Aharonov-Casher effect&HL and the Sagnac phase 
shift 11 - have been studied. Although many of these fea- 
tures can be understood in terms of a single particle pic- 
ture, quantum rings are also attractive laboratories, be- 
ing quasi-lD systems, to study many-body effects. A 
review about this topic can be found in Ref. 1 12L In con- 
junction with spin-orbit interaction, quantum rings are 
an active field of research for spin manipulation 13 . Pro- 
posed devices, based on a ring geometry, include spin- 
filtersii, spin beam splitters^, spin current parametric 
pumping^, spin qubitsi 7 -, or quantum rings which can be 
used to compose networks for spintronic applications^. 

Topological insulators (TI) are a new class of time- 
reversal invariant materials that have a bulk gap and pro- 
tected surface states. These materials are distinguished 
from a trivial insulator (with no topological protected 
surface states) by a Z2 topological number—. In two 
dimensions, these materials exhibit the quantum spin 
Hall effec t 20 ' 21 and typically are small gap semiconduc- 
tors, with the peculiarity that their gap, originated from 
strong spin-orbit coupling, is inverted with respect to 
that of a normal insulator, allowing for the presence of 
topologically protected edge states22r— . At each bound- 
ary two gapless counterpropagating edge modes exist 
that have opposite spin and are related by Kramer's the- 
orem, known as helical edge modes. A successful realiza- 
tion of a two dimensional TI phase has been obtained in 
HgTe/CdTe quantum wells, where the topological phase 
has been confirmed with the observation of the ballistic 



helical edge channel a 23 ' 24 . TIs, with their dissipation- 
less edge states, represent a new promising route to spin 
manipulation in semiconductors. 

A ring geometry in a TI with its well localized wave- 
functions is optimal for spin manipulation, using a mag- 
netic flux. It has been proposed, on the basis of numerical 
transport calculations, that the AB effect in a topological 
insulator disc can be exploited to perform spin rotation 
and spin filtering 2 ^, or to create a topological spin transis- 
tor—. Similar properties are expected in a general closed 
disk geometry connecting two quantum point contacts in 
which coherent oscillations in the magnetoconductance 
are expected due to the AB effect 2 -^. It is clear that, if 
experimentally confirmed, these properties can make the 
disk, ring and hole geometry extremely interesting also 
from an application point of view. 

In this paper we consider the boundstates of a ring 
made from a two-dimensional TI (see Fig. [I}. The bound- 
states are found analytically as a function of a magnetic 
flux, piercing the ring hole, and dimensions of the ring in 
the topological non-trivial (edges states and bulk states) 
and topological trivial (only bulk states) regime. We 
show that the specific bandstructure of such TIs allows 
for manipulations of electron spin states using the mag- 
netic flux and their read-out through a spin-selective per- 
sistent current in the ring. 

In addition, the special topology of a ring introduces 
a new type of two-level system consisting of inner and 
outer-edge localized bound states where their mixing can 
be tuned by the flux. Together with the spin of an elec- 
tron, these features allow for new ways to coherently ma- 
nipulate quantum states in topological insulators. 

The paper is organized as follows. In section II we 
solve the edge and bulk states of a TI in cylindrical co- 
ordinates, described by an effective four-band model for 
HgTe/CdTe quantum wells 21 , exactly. We summarize 
here, for the benefit of those who are not intersted in the 
more mathematical details, the main results of section II. 

The total out-of-plane angular momentum is a con- 
served quantity of the system and we associate the half- 
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FIG. 1: (Color online) Scheme of a topological insulator ring 
of radius a and width w. A magnetic field, generating a flux 
$ is threading the ring. 



odd integer quantum number m to it, this fact permits us 
to separate azimuthal and radial variables in the spinor 
[Eq. ©]. We can express the radial part of the spinor 
components [Eq. (|22p and (|2T>|) ] in general as a combina- 
tion of ordinary and modified Bessel functions, both for 
the bulk [Eq. (gj)] and edge modes [Eq. (jgTJj) ]. 

We find the eigenstate spectrum by assuming vanishing 
boundary conditions for a ring, a disk or a hole geometry, 
obtainining transcendental secular equations, that we nu- 
merically solve for the eigenenergies, and eventually we 
obtain the corresponding eigenfunctions. So in general 
an eigenstate of the system can be labelled as |n,m,r), 
with n the radial quantum number, associated with the 
number of nodes of the radial wavefunction, the total 
angular momentum m and the spin r = ±1. 

In Section III we show the results of our calculation, in 
particular focusing on the bound states of a ring subject 
to the AB effect. We present the edge dispersion for a 
disk, a hole and a ring, where a minigap arises due to the 
confinement, which exponentially shrinks with increasing 
ring thicknesses. We show that the sign of the Dirac 
mass, which distinguishes the topologically trivial from 
the topologically non-trivial insulator, is important to the 
eigenstates in the gap but also in the bulk spectral region. 
We also calculate the persistent currents as a function of 
the magnetic flux. 

In Section IV we consider the inclusion of Rashba spin- 
orbit interaction which couples the two spin-blocks. The 
main effect, for a sufficiently small interaction, is the 
opening of a spin splitting gap at specific values of the 
magnetic flux. As a consequence, we can perform coher- 
ent spin manipulations with the magnetic flux. 

In Section V we present a scheme for the coherent ma- 
nipulation of quantum states stored locally in the spin or 
the edge degrees of freedom of a TI ring, by means of an 
AB flux. 



II. FOUR BAND MODEL OF HG-TE 
QUANTUM WELLS 

For a realistic description of a two-dimensional 
topological insulator, we use the effective four-band 
model 21,24 for HgTe quantum wells, composed out of the 
Ei and H\ subbands, and the electron spin degree of free- 
dom |±) (see Appendix IB"! for a discussion of these basis 
states and the meaning of spin). The four-band Hamil- 
tonian is built with respect to the basis {|Ea+), |-ffi+), 
\Ei-), \H\— )} (see Eq. (jBll) ). where the E\ states are a 
mixing of the s-like Tq band with the Tg light-hole band, 
and Hi represents basically the Tg heavy-hole band. In 
this basis, ± stands for the sign of the total angular mo- 
mentum projection of the E\ and H\ bands, i.e. the sum 
of spin and orbital angular momentum. In the absence 
of inversion symmetry breaking terms (as the Rashba in- 
teraction discussed in Section IV), the Hamiltonian is 
separable into two blocks related by time-reversal sym- 
metry 21 

g_(h(k) \ 

a ~ \ h*(-k) J w 

with h(k) = e(k)l 2X 2 + d a (k)a a and d a {k) = 
(Ak x , —Ak y ,M(k)) where a a is the vector of Pauli ma- 
trices. In Eq. ([T]), e(fc) = C — Dk 2 , and the Dirac mass 
term M (k) = M - Bk 2 with k 2 = fc 2 + fc 2 . 

The parameters A, B, C, D, M depend on the QW ge- 
ometry. Realistic estimates can be made by a comparison 
of the effective model with a well established 8x8 Kane 
Hamiltonian 28 . In the following we will use the param- 
eter values A = 375 meV nm, B = —1.120 eV nm 2 and 
D = —730 meV nm 2 . The parameter C results only in a 
shift of the energy and we can set C = without loss of 
generality. The Dirac rest mass M depends on the QW 
thickness and M < corresponds to the inverted (topo- 
logically non-trivial) regime whereas M > corresponds 
to the normal (topologically trivial) regime. 

This model is able to describe both the edge states and 
the bulk states in a reasonable energy ranged. We note 
that Eq. (HJ) is not restricted to HgTe QWs but describes 
also other two-dimensional TIs like thin films of Bi2Te3 
and Bi2Se3. 30 For the ring, it is useful to express Eq. ([1]) 
in polar coordinates x = r cos 9, y = r sin 9. 

We now use k — — jV and perform the unitary trans- 
formation 

U = i [(r - T z )a y + (r + T z )a ] , (2) 

which changes the old basis into \Hi+), 
—i\Hi—), i\E\— )}, we obtain the transformed Hamilto- 
nian, 

H = C + Mr z a z + (Dl + Br z a z )A - iAe irT * e Tl (3) 
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U(r,t 
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r 2 dO 2 



or r oO s 



(4) 
(5) 



where r a , <r a represent the Pauli matrices for the spin 
(Kramer's partners ±) and pseudo-spin (El, HI) com- 
ponents, respectively. 

In the following table we define the orbital angular 
momentum l z , the pseudo spin operator S z and the time 
reversal operator T for the unsymmetrized and the sym- 
metrized basis. 



H 


H 


l z = -ihdg 

S z — 2 Tz< -^ z 

T = ITyCT^K 


l z = -ihde 
S z — ^T a z 

T = -T X <7 y K 



(6) 



We further introduce a total angular momentum oper- 
ator as j z = l z — S z which commutes with H 



[H,J Z 



0. 



(7) 



The eigenstates of j z therefore can be used for the ring en- 
ergy eigenstates. Consequently, we can write the spinors 
for the two blocks (r = ±1) of H as 



# m>T (r,0) = 



Am9 



/2vr 



m.T / 

Xi 

m.T/ \ - 

X2 ( r )e 



where it is easy to verify that 



hm^>r. 



(8) 



(9) 



Now we need to find the radial wave functions Xi 2( r )- 
We note that the diagonal part of H (putting II = 0) is 
solved by Bessel functions. Similarly, if we put A = 0, 
the radial solution for each of the two blocks of H can 
be expressed in terms of Bessel functions^. We con- 
sequently use the ansatz xT' T ( r ) = /m+±(-^ r )i with 
/ proportional to a function of the Bessel family, de- 
fined below (this will result in x™' T (r) oc f m _i(Kr) 
as we will show below). To include a magnetic flux $ 
threading the ring hole, we consider the vector potential 
A = ($/27rr) eg. Using the minimal coupling procedure 
—ihW — > —ihV + eA, with — e the electron charge, the 
inclusion of the magnetic flux is equivalent to 



-> de + i — 
$ 



(10) 



with the magnetic flux quantum $ = h/ e - The radial 
wave functions satisfy the Bessel differential equation 



A™ ± i/™±i(ifr) = -K 2 f,- n± i(Kr), (11) 



with 



d 2 Id 1 /_ 1 
or £ r or r z \ 2 



(12) 



where m = m + $/$o- Using H^ m T — E^> m ^ Tl we arrive 
at the following equations 



.R 
K 



d r 



Xl =X2 



(13) 



(E+DK 2 ) 2 -(BK 2 -Mf 

T 2 



with 



R 



AK 



(BK 2 - M)t -(E+ DK 2 ) 



xT T = o, (w) 



(15) 



Eq. (|14j) is satisfied, employing the assumption of xT 
being a Bessel function (see Eq. (fTTj) ). if 



A Z K A = (E + DK 2 ) 2 - (BK' Z - M) 



(16) 



Therefore, for a given energy E, there are four possible 
values for the wave vector K 



K 2 



-F±^F 2 -Q 2 , 



with 



F 



A 2 - 2(BM + DE) 

2{B 2 ~- 
M 2 - E 2 



D 2 ) 



B 2 -D 2 



(17) 

(18) 
(19) 



Solutions with K 2 > correspond to propagating 
modes, and f(r) is therefore described by the Hankcl 
functions Hu and For K 2 < 0, instead, / can be 

expressed with the modified Bessel functions K v (\K\r) 
and /j/d-ftTlr) describing an evanescent behavior (expo- 
nentially growing and decaying) along the radial direc- 
tion. 



A. Condition for bulk-propagating modes 

In the spectral range inside the bulk bands (\E\ > 
\M\), we expect Eq. (fl?]) to produce two real (K = ±Ki) 
and two imaginary K = ±iK 2 solutions. Note that for 
F 2 — Q 2 < 0, also complex solutions for K 2 are possible. 
However, F 2 — Q 2 > as long as M > ^- which is the 
case for HgTe QWs. Here K 1 (K 2 ) are the modulus of the 
real(imaginary) solutions of K obtained with the +(— ) 
sign from Eq. ([T7|) . The required condition is that 



F < y/ F 2 - Q 2 



which is satisfied if 



Q < 0. 



(20) 



(21) 



In HgTe QWs, \B\ > \D\, and the condition is automat- 
ically satisfied if \E\ > \M\, as expected. 
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Within this spectral region, the first component of the 
spinor in Eq. (JSJ can be expressed as 



X?' T (r) = hF m+ i(r), 



(22) 



with 



h = (h 1 ,h 2 ,h 3 ,h 4: ) (23) 
F v {r) = fi?W(ei),^ 2) (6),^(6),^(6)), (24) 



where £j = JQr, i = 1,2. The expression for the second 
component can be found from Eq. (|13[) and using the 
following property of the Bessel functions 



d r g v {Z)±-g v (Z) = ±riKg v - 1 (Z), 



(25) 



with rj = +1 for Hu \ and I„, and r\ = — 1 for if„ 
We obtain 



(26) 



x ™> T (r)=ihRF^(r), 
with the diagonal matrix R given by 

R = diag (rQ,rQ,R%1-r£! 
where, according to Eq. (|15p . we define 



R 



(t) 



(±M t 2 - M)t -(E± DKf) 



(27) 



(28) 



with ± depending on whether K i is the modulus of a real 
or an imaginary solution of Eq. (|17p . 



B. Condition for edge modes 

In order to have four imaginary solutions of Eq. (fTT| , 
in the region \E\ < \M\, we need, being \B\ > \D\ and 
Q 2 > 0, to satisfy the conditions F > and F 2 - Q 2 > 0. 
The latter turns out to be the more restrictive one. In 
particular, if we want to observe the Dirac point of the 
edge states (of typical value E = —A1D/B), the material 
parameters have to satify the following inequality 32 



D 



A 2 AM 

^Td 2 > IT" 



(29) 



Again, Ki and K 2 are the wave vector moduli of the 
four modes with imaginary wave vectors obtained from 
Eq. ([17} with the + and — sign, respectively. 

The spinors for edge states are also obtained from 
Eq. © with components described by Eqs. (|2"2"|) and (|21)1) . 
but with an alternative definition of the basis vector F v 
and the diagonal matrix K. They are given by 



F v {r) = (^(a),^(6),^(6),^(6)) 



(30) 



diagfiZ^ -4^,^.-4:1). (31) 




FIG. 2: (Color online) (a) Upper (U) and lower (L) bands 
of the edge states for a ring with inverted Dirac mass M = 
— 10 meV, radius a = 500 nm and width w = 200 nm for spin 
up and down blocks with solid (r = 1) and dashed (r = —1) 
lines. Squares and circles represent the internal and external 
boundary (helical) edge states calculated independently, filled 
for t = 1 and empty for t = — 1, which would result from a 
very large width w where no minigap develops. They could 
also describe a hole and disk geometry, (b) Radial probability 
density p(r) = I^MI 2 of the edge states of the upper band 
with t — 1, for fh = and fh = ±|. 



C. Boundary conditions and secular equation 

In order to obtain the energy levels and bound states 
of the ring, we need to solve the boundary problem. We 
assume a ring of radius a and width w and impose van- 
ishing boundary conditions at r = a ± | 



(32) 
(33) 



Details of this calculation can be found in Appendix [X] 
Here we just note that the four former conditions plus 
the normalization of the wave function fix h and lead 
to a transcendental secular equation for E. We here re- 
mark in connection with rings in graphene 3 ^, that the 
inclusion of quadratic terms allows for a formulation of 
vanishing boundary conditions which is not allowed in 
the Dirac equation (only linear terms in k) of graphene 
as the wave function would vanish identically. There, a 
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FIG. 3: (Color online) Minimal and maximal values of the 
minigap (it is a periodic function of the magnetic flux) for 
a ring with radius a — 500 nm as a function of w between 
100 nm and 300 nm. 
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FIG. 4: (Color online) Low-lying bound states corresponding 
to m = 5, t = 1, a Dirac rest mass of M = ±10 meV, ring 
radius a = 500 nm and width w = 100 nm. The change of 
sign of the mass term is accompanied by the disappearence of 
the edge states and by a non-trivial shift of all eigenenergies. 

confinement on a circle can be formulated with an infinite 
mass boundary conditio n 31 i 33 . 

For the bulk-propagating region we obtain the secular 
equation described by Eq. (|A1|) , while for the edge states 
we obtain Eq. (|A2I) . For a fixed m and r, the secular 
equation is numerically solved. In the bulk-propagating 
region we obtain the eigenvalues E n _ m ,T where n = ±2, 
±3, . . . represents a radial quantum number, positive for 
conduction band modes and negative for valence band 
modes. The corresponding eigenstates are given by the 
spinors ^ n ,m,T or in ket notation \n, m, r). We can think 
of the edge states, when system parameters allow their 
existence, as the system's eigenstates with n = ±1. 

III. BOUND STATES AND PERSISTENT 
CURRENT OF THE RING 

When the mass parameter M is negative (the topo- 
logical non-trivial regime), our model accepts solutions 
in the (inverted) gap region E £ (M, —M). This solu- 
tions are exponentially localized near the external and 
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FIG. 5: (Color online) Dispersion curves for the first bulk- 
conducting mode n = 1, as a function of the magnetic flux for 
a ring in the inverted regime (M = — 10 meV) of radius a = 
100 nm and width w = 75 nm (note that this mode originates 
from the edge states pushed into the bulk conduction band 
by their overlap). States for m = ±1/2, ±3/2 and ±5/2 are 
shown, in (a) for the block r = 1 only, and in (b) for both 
blocks r = ±1 with continuous and dashed lines, respectively. 



internal boundaries r — a ± ^ respectively, and are ad- 
dressed as edge states. In Fig. [^a) we show the edge 
state dispersion for a relatively large and wide ring of 
radius a = 500 nm and width w = 200 nm calculated 
with Eq. (|A2I) for a boundary condition including both 
ring edges. The physical content of this picture can be 
undestood remembering that for a given magnetic flux 
the eigenstates of the ring are found at fh = m ± $/$o 
with m an half-odd integer. It means that with $ we 
can move the eigenstates along the dispersion curves in 
Fig. Uta), always having two subsequent states at uni- 
tarily spaced intervals on the m-axis, i.e. for $ = we 
have states at m = ±1/2. The curves with squares and 
circles represent the internal and external edge states re- 
spectively, as calculated from the decoupled boundary 
conditions with Eqs. (|A8|) and (IA5[) . We can also think 
of them as originated from edge states of a disk (external 
one) and of a hole in a plane (internal one). Their dis- 
persions are different because they are centered on two 
different effective radii (approximately at r = a ± |). 
We note that their slope (which is an energy) is in fact 

approximately given 25 by ^ J B2 g-F 2 '. This emerges nat- 
urally from the fact that the radial wave function has to 
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FIG. 6: (Color online) Spin-dependent persistent currents cal- 
culated for N = 1, 2, 3 and 4 electrons in the ring's first 
bulk-propagating mode as a function of the magnetic flux for 
a ring of radius a = 100 nm and width w = 75 nm. In (a) 
only the spin block r = 1 is occupied, while in (b) electrons 
can occupy both spin-blocks. 



close on itself in a roundtrip along the circumference. 

When the ring width is small enough to permit a fi- 
nite overlap of the helical edge states, internal and ex- 
ternal edge states get mixed. Consequently, as shown in 
Fig. rjja), the dispersion curve is modified into the anti- 
crossing upper and lower bands separated by a minigap, 
as it was also predicted for a straight TI ribbon of finite 
width 32 . In Fig.rjfb), we show the radial probability den- 
sity of the edge state spinors of states at different points 
along the upper band (U) for r = 1. This picture gives 
us an understanding of the localization properties of the 
edge states near the minigap. The state at fh — (at the 
anticrossing) is equally shared between the internal and 
external edges. Away from this anticrossing point (we 
plot the states with fh = ±|), where the dispersion curve 
resembles that of the uncoupled edges, the edge states are 
strongly localized near one or the other boundary. This 
means that we are able to drag an electron from one edge 
of the ring to the other by tuning $. Conversely, for the 
lower (L) band fh = ±^ states exhibit opposite local- 
ization properties as compared to the upper band. The 
implications of this tunability for quantum information 
processing will be discussed in Section V. 

For certain values of the flux (when a state is in corre- 
spondence with fh — 0), as expected, the minigap expo- 
nentially shrinks with increasing w (see Fig. [3]). However 



in a fully confined system, as our ring, the discreteness of 
levels can also lead to a saturation of the minigap at a fi- 
nite value. For sufficiently small w the edge modes can be 
pushed into the bulk-region, losing their edge-character 
and becoming the first of the bulk-propagating modes. 
In this case, the minigap is comparable or larger than 
the bulk gap itself (i.e. larger than 2|M|). 

The presence of the edge states with energy in the bulk 
gap depends on the negativity of the Dirac rest mass 
M, while for positive M only bulk states exist. But the 
change of the mass sign has also effects on the bulk en- 
ergy levels as can be seen in Fig. 21 which shows the 
lowest lying eigenergies with fh — ^ and r = 1 for a 
ring of radius a = 500 nm, width w = 100 nm and for 
M = ±10 meV. Therefore the presence of the topological 
insulator phase in confined systems can also be deduced 
from bulk-propagating properties as measured by optical 
spectroscopy or tunneling transport experiments through 
the ring. 

In Fig. [51 we show the dispersion curve as a function of 
the magnetic flux for the first bulk-conducting mode for 
a relatively small and narrow ring of radius a — 100 nm 
and width w = 75 nm. We select a ring of small radius in 
order to enhance the magnitude of the persistent current, 
which is proportional to the derivative of the dispersion 
curves with respect to the magnetic flux, which grows 
with decreasing a. The small confinement length w leads, 
instead, to well separated eigenvalues corresponding to 
different radial quantum numbers n. The energy level 
spacing, between energy eigenvalues of different radial 
quantum number n (see Fig [5]), scales, in fact, with 1/w. 
In Fig. [5ja) we show the states of the spin block r = 
1 with m = ±1/2, ±3/2 and ±5/2. As can be seen, 
also in the absence of the magnetic field, effective time 
reversal symmetry (m — > —m) is broken within a single 
block. The presence of the second spin-block restores 
time reversal symmetry when $ = 0, and states with 
the same n, and opposite r and m are time reversal (or 
Kramer) partners. The dispersion curve is periodic in $ 
with periodicity given by the flux quantum $o- 

The persistent current at zero temperature is given by 



E 



dE n , 



<9$ 



(34) 



where the sum runs over all occcupied states. In Fig. [5] 
we show the persistent current for a ring with the same 
parameters as used for Fig. [5j We consider an electron 
occupation number TV up to 4 occupying the first radial 
mode with positive energy, subtracting the current con- 
tribution coming from all states below. In Fig 13 a) we 
consider injection of spin-polarized electrons correspond- 
ing to r = 1. As a consequence, the persistent current is 
finite even in the absence of a magnetic flux, and its sign 
directly reflects the spin polarization. While in Fig^b), 
both blocks can accommodate electrons and, accordingly, 
the current vanishes at zero magnetic field, respecting 
the time reversal symmetry of the model. However, we 
note that the spin current can be finite at $ = 0. The 
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FIG. 7: (Color online) Details of the dispersion curve of the 
first radial mode in the presence of the Rashba term (qo = 
20 meV nm) for a ring of radius a = 100 nm and width 
it) = 75 nm. Lines are dispersion curves without the Rashba 
term and points after its inclusion (a thin line is used to link 
the points to emphasize the spin-splitted bands). The Rashba 
term induces a splitting at the crossings of with 
1, —1), and of |1, — i 1) with |1, |, —1) (not visible in this 
range of $). 



currents are periodic in $ with the same periodicity as 
the dispersion curves in Fig. [5j and show the presence 
of kinks whenever different bands cross and electrons are 
forced to move from one to the other in order to maintain 
the groundstate of the system. 

The persistent currents in Fig. [B] are calculated for an 
edge state (n — 1) pushed into the bulk conduction band 
by the overlap between internal and external edge states. 
But generally, a similar behavior is observed for other 
bulk states of the ring (\n\ > 1, not shown), where a 
spin-selective behavior is not as obvious as in the contri- 
bution from the helical edge states. The reason for this 
phenomenon is that the Dirac mass M(k) breaks the sym- 
plectic symmetry p — > —p and a — > —a in a single block. 
In graphene rings with a smooth ring boundary the same 
symmetry is broken and leads to a valley- dependent per- 
sistent current^. This physics in principle would allow 
to use rings in topological insulators as spin-detectors 
where the direction of spin is mapped to the direction of 
the persistent current. Indeed currents shown in Fig^a) 
and (b) are well in the range accessible to experiments 8 . 
The persistent current has also been calculated for a disk 
in an homogeneous magnetic field 34 . However, the topol- 
ogy of that system is qualitatively different from a ring. 



IV. RASHBA SPIN-ORBIT INTERACTION 

We consider the inclusion of a Rashba spin-orbit inter- 
action as derived in Ref. l35l . retaining only terms to first 
order in k, which gives the most dominant contribution 

Hr = « + 2 ~ ( T y kx ~ Txk y^ ■ ( 35 ) 



We use cylindrical coordinates and transform it into the 
symmetrized basis using Eq. ([2]). We obtain 

Hr = ta — — e lT * a I Tyd r -T x — I , (36) 

where ao is the Rashba coefficient. In particular, the 
Rashba interaction in HgTe quantum well can be tuned 
with a gate voltage, so that the induced spin-splitting 
gap varies from zero up to tens of meV— . 

We note that Hr does not commute with H and j z , 
and therefore will in general mix energy eigenstates of 
H . We can still find a conserved physical observable by 
noting that the only effect of Hr is to couple \E\, +) with 
\Ei,—). This motivates us to define the total angular 
momentum restricted to the i?i-band 37 

^0=^)+^ (37) 

where the operators are the same as described in Ta- 
ble ([HI), but their action is now restricted to states of 
the Ei band only. This can be obtained by applying the 
projector _p( £l ) = CT "+J* cr * to the operators defined in Ta- 
ble so that H El) = i z P^ and sf l] = S Z P^\ It is 
clear that the Rashba term Hr must conserve ji El ^ and 

A (E ) 

indeed [HR,j z 1 ] = 0. On the other hand the spinors 

^n.m,T are eigenstates of H and eigenstates of jz with 
eigenvalue m + r 

] z El) ^n. m ,r = h(m + T) V n , m ,T. (38) 

Therefore, Hr can only couple eigenstates of H from the 
blocks t = 1 and r = — 1 with 

m (T= _i) - m (r=1 ) = 2. (39) 

The conservation of the "i?i-band" total angular mo- 
mentum imposes that the only non-trivial matrix ele- 
ments exist between the spinors \n,m, 1) and \n',m + 
2, —1), with n and n' two quantum numbers of the radial 
quantization, and m the quantum number of the total 
angular momentum. In the following we will assume the 
eigenvalues of two different radial eigenstates (n ^ n') 
to be sufficiently distant in energy, and the Rashba term 
sufficiently small, such as to make the mixing induced by 
spin-orbit negligible. Actually the strength of the Rashba 
term can be tuned in topological insulators, while the dis- 
tance between neighbouring radial eigenstates increases 
for smaller w, so that we can always restrict ourselves in 
a parameter range which makes the previous assumption 
valid. 

The relevant matrix element of Hr is 

H R ( n , m ) = (n,m + 2, -l\H R \n, m, 1) = 

f, (m+2,-1)* f a (m, + l) , AM 

= ia / drrx\ \o r J x\ , (40) 

which, using the Bessel function's property Eq. (|25p . can 
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FIG. 8: (Color online) Schematic view of a flux tunable two- 
level system. |0) and |1) states refer to internal and external 
edge states, when the splitting huj is the minigap. Another 
definition is |1) = |r = +1) and |0) = |r = —1), when we 
are dealing with the spin degree of freedom and the energy 
splitting is due to the Rashba interaction. 



be put in the form 



H 



R(n,m) 

tag / drr 



IMF, 



|n,m+2,-l) 



hKR 



rn+4 



|n,m,l) 



(41) 



with K = diag (Ki,r]Ki, K2, —K2), and where the coef- 
ficients h, R and K, in squared parenthesis, are calcu- 
lated for the indicated spinors |n,m,T). This expression 
is valid for both bulk and edge states, for which r\ = ±1, 
respectively. 

In Fig. we show the effect of a Rashba term with 
ctQ = 20 meVxnm on the dispersion relation of the first 
bulk modes of a ring with the same parameters as in 
Fig. [3] We note that the most significant signature is 
the opening of a splitting with an anticrossing behav- 
ior between the states |1,— and |1, |,— 1), and be- 
tween their time reversed ones |1,— 1,1) and |1,|,— 1) 
(not shown). 



V. COHERENT MANIPULATION OF SPIN 
AND EDGE DEGREES OF FREEDOM 

An AB-ring in a two-dimensional topological insulator 
can be used to process quantum information. Our idea is 
to use the magnetic flux to manipulate a two-level system 
of states |1) and |0), which are good system eigenstates, 
but which can be coupled by some interaction for spe- 
cific values of the flux (see the scheme in Fig. [HJ . At res- 
onance, corresponding to the operational magnetic flux 
<f> c , the eigenstates of the system can be described by 
\U) = (|l) + e*|Q))/V2 and \L) = (|1) - e*|0)) /y/2, 
for the upper (U) and lower (L) band. 

We identify in our TI ring two pairs of such two-level 
systems. The first one, which is specific to a topological 
insulator ring, makes use of the external ( 1 1) ) and internal 
(|0)) edge states, which realize a system similar to that 



of a charge-qubit in a double quantum dot^S, but where 
charge separation is due to the nature of the edge states. 
A ring has actually an even number of Kramer's part- 
ner and therefore it is topologically trivial, in the sense 
that scattering is allowed between the two different set 
of helical edge states localized at the internal and at the 
external boundaries. However for sufficiently large w a 
degree of the topological protection remains, because the 
scattering is suppressed by their localization properties. 
As shown in Fig.[^b), at zero flux, states of the same spin 
r with m = 1/2 and —1/2 are strongly (exponentially) 
localized at the external and internal boundary respec- 
tively. Let us imagine to inject a state of definite spin r 
into the state m = 1/2 which is strongly localized at a 
single edge (see Fig. [5]) and then non-adiabatically vary 
the flux from zero to the operational value $ c = — $o/2, 
corresponding to fh — 0. In this situation the eigenstates 
of the system are equally shared between the two edges 
and an energy splitting huj (the minigap) arises, due to 
the inter-edge overlap. The original state will now be 
subjected to Rabi oscillations between the internal and 
the external boundaries with angular velocity uj. Read 
out of this two-level system could be performed by con- 
necting separately, with tunneling contacts, the internal 
and external boundary, or measuring the sign of the flux 
generated by the persistent current (cf. Fig. [5]). 

Another, more canonical choice, initially proposed in 
Ref. HI is to use the spin degree of freedom and to asso- 
ciate |1) = |+) and |0) = |— ) with ± referring to the up- 
per and lower block of the TI Hamiltonian with r z = ±1. 
This two-level system is indeed topologically protected in 
a ring for a sufficiently large w, or in a disk or hole geom- 
etry, from time-reversal symmetric scattering, but still is 
subject to external decoherence sources, like the hyper- 
fine interaction with an underlying nuclear spin system. 
The tunable Rashba term can offer a way to manipu- 
late the spin-degree of freedom^!!, as implied by Fig. [7J 
where a spin splitting gap (tko) opens between the state 
|+) = I m, t z = +1) and the state |— ) = \m + 2, t z = —1) 
around $ = $0/2. The spin state could be observed by 
detected the persistent current, which is opposite for the 
two spin states (cf. Fig. 

Both spin and edge (i.e. charge) two-level systems 
can be operated by the magnetic flux. For example, 
let us imagine to be initially in state |1), in a range of 
the flux value for which it is a good eigenstate of the 
system (see Fig. [8]). Swiftly turning the magnetic flux 
to the operational value 4> c , and keeping it there for a 
time t, will lead to the final rotated state described by 
cos(wt/2)|l) +ie-^sin(w£/2)|0). We can switch off the 
rotation by restoring the magnetic flux to its initial value. 
For the stability of the two-level system it is required that 
the dephasing and relaxation time are much larger than 
the Rabi period 2tt/uj. For these two kinds of two- level 
systems (edge and spin) , the splitting can be on the order 
of 1 meV for reasonable parameters a and w, correspond- 
ing to a Rabi period of 4 ps. 

One can also exploit the ability to control the localiza- 
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tion properties of an edge electron, that is, the possibil- 
ity to move it from the internal to the external bound- 
ary by means of the flux, in a ring network. In such a 
system, spin qubits can be stored in the edge states lo- 
calized on the internal boundary, the position in which 
they are protected from the interaction with the network. 
Bringing them to the external boundary, corresponds in- 
stead to put them in an active position, allowing for their 
measurements or for the interaction with qubits on other 
rings. 

VI. CONCLUSION 

In the present paper we have solved the eigenstates of 
a topological insulator for a ring, a disk and a hole struc- 
ture, taking HgTe quantum wells as a practical example. 
We solved analytically the four-band model with both 
linear and quadratic terms, yielding a transcendental sec- 
ular equation which we solved numerically as a function 
of the magnetic flux threading the ring. We calculated 
persistent currents, which can be exploited to measure 
the spin state of the system. The effect of a Rashba 
spin-orbit term has also been taken into account, show- 
ing that its main effect is the opening of a spin-splitting 
gap at specific values of the magnetic flux. We proposed 
that, in a TI ring, spin and edge (charge) degrees of free- 
dom, where either up/down spins or internal/external 
edges are used as quantum basis states (two-level sys- 
tems), can be coherently processed. Rashba interaction 



can be used to manipulate spin, while a finite edge-state 
overlap can be exploited to manipulate the edge degree 
of freedom, whereas the AB-flux is used to introduce co- 
herent dynamics of the two-level system. 
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Appendix A: Boundary conditions and secular 
equations 

1. Bulk-propagating modes 

The bulk-propagating eigenstates of a TI ring, of total 
angular momentum m, can be written as as a linear com- 
bination of the spinors of the four modes K = ±K\ , and 
±iK 2 as described by Eqs. ©, (JUJ) and ([21]). Now, van- 
ishing boundary conditions (Eq. (1321) and (1331) ) impose 
the spinor components to vanish on the internal and ex- 
ternal border and, toghether with the normalization con- 
dition, fix both the four coefficients h n and determine a 
the secular equation for the system eigenvalues. We ob- 
tain the following secular equation 



z 



= 



R 2 A_(K + (H 2 , 2 )H + (£ 1A ) - K+fedH+faj)) - i? 1 A + (^_(6, 2 ) J ff-(6,i) - #_(6,i)ir_(£i >2 )) 



(Al) 



where we adopted the following notations 





= K u r v ; 




U, V 


= 1,2 




w 

a ~ 2"' 




r 2 


w 

= a+ 2 


9± 


= Sm±§; 




9 = 


H, K,I 


Ai 


= 


)K- 




+ /_(6a)iM6,i) 


A 2 




)K- 


(6,2) 


+ I_(&, a )JT + (£ a , 2 ) 


A+ 


= 


)K, 


-(6,1) 


-/+(&,!)#+ (k,2) 


A- 


= I- (6,2 


)K- 




-/_(6,l)^-(6,2) 



Once the transcendental equation Eq. (jAII) is solved for 
the eigenenergy, the coefficients h are obtained by solving 
the boundary conditions. Note that these results depend 
on r through the terms R\ and R 2 . 

For a given value of the quantum numbers m and r, the 
secular equation gives a series of eigenstates labeled by a 



radial quantum number n, with \n\ corresponding to the 
number of nodes of the second component \ 2 . Generally 
with M < edge states have n = ±1 and bulk propagat- 
ing states are labeled by n = ±2, ±3,. . . , where ± are for 
conduction and valence band states, respectively. When 
M > no edge states are present and bulk-propagating 
modes are labelled by n = ±1,±2,. ... As an example, in 
Fig. [9l we show the numerical solution for rh = of a TI 
ring in the inverterd regime (M = — 10 mcV) of radius 
a = 1000 nm and width w = 50, 100 and 200 nm. The 
red curve corresponds to $s(Z) in Eq. (|A1[) and its zeros 
define the system eigenvalues. Note that the first zero, 
for the case with w = 50 nm, corresponds to an edge state 
pushed into the bulk conduction band by the overlap be- 
tween internal and external edges (and it is therefore the 
n = 1 state). The dotted lines correspond to the an- 
alytic fit sin [Ki (w - K2 1 )] where K h2 are calculated 
at the energy E + M^, with the parameter t « 50 nm. 
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E (meV) E (meV) E (meV) 



FIG. 9: (Color online) Q(z) defined in Eq. [fAT]). whose zeros give the eigenvalues of a ring with radius 1000 nm and width 50 
(a), 100 (b) and 200 nm (c). The dashed curve is an analytic fitting curve, described in the text. 



As can be seen, this fitting curve well reproduces the 
spacing between the eigenvalues of the system. The solu- 
tions are therefore determined by the energy E for which 
Ki(w — K 2 ) = mr. With different values of the mass, 
we need to adequately fit a new £, which appears to de- 
crease with M. For M = 5 and 15 meV the best fit is 
given by I = 90 and I — 37 nm, respectively. 

In Fig. [TUta) and (b) we show the n = 1 and n = 3 
radial modes of the system described in Fig. EJa), cor- 
responding to an energy of 17.2 meV and 96.8 meV, re- 
spectively. 



2. Edge states 



In the gap spectral regions of a TI ring, we find evanes- 
cent solutions (i.e. helical edge states) if M < 0. For a 
set of quantum numbers m and r, Eqs. ([5]), (|2"2"|) and (|26p . 
togheter with the definition of Eq. (j30]) , describe the two 
spinor components. Now, solving vanishing boundary 
conditions (Eq. (|32l) and (|3"3")l ) and the normalization of 
the wavcfunction, we are led to the secular equation 



Q 
Q 

z 



(A2) 



A 2 + R 1 K+^ 1 , 1 )I-^ 1 , 1 )) - A 1 (M_(£ 2 , 2 J + (a, 2 ) + M+(a,2)I-(6,2)) 



r 2 a„ (*r + (k, 2 )/+(6,i) 

A 2 (R 2 K + (Z ltl )K-(Z 2>1 ) - 



- X+(6,i)/+(Ci,2)) - RiA+ (A-_(€a, 2 )/-(Ci,i) - tf-(£a,i)M&,2)) 



i? 2 A_ {K + (£ ltl )K + farf - K+farfK+fa,!)) + RiA+ (K- (£i,i)K- (6,2) - ^-(6, 2 )^-(6,i)) 



This equation fixes the eigenenergies, while the corre- 
sponding spinor coefficients h can be obtained by solving, 
for a fixed energy, the boundary conditions. 

a. Isolated boundaries: disk and hole 

If we consider w sufficiently large so as to effectively 
isolate the two boundaries, or just consider a disk or hole 
geometry, we can independently solve the internal (— sign 
in Eqs. ([22]) and ([22])) and external (+ sign in Eqs. ([3"2]) 
and (|33p ) boundary conditions. The internal boundary 
conditions (hole geometry), for themselves, are compati- 
ble with h 2 = /14 = and give 

/+(6,a) + h 3 I+(&,2) = (A3) 
Ril-(&,i) + h 3 R 2 I-(&,2) = 0, (A4) 

which results in the secular equation 

/+(a, 2 ) _ Ri 1-^1,2) rA ,^ 

M&, 2 ) R 2i- (6,2)' 1 J 



For a given energy, the coefficients hi and /13 are then ob- 
tained with the same boundary conditions. The problem 
of holes in a TI has also been addressed in Ref. |4lJ. 

For the external boundary (disk geometry) , consistent 
with h\ = h% = 0, we find 

^+(6,i) + Mv + (£ 2 ,i) = (A6) 
flx^feO + ^M-fei) = 0, (A7) 
and, therefore, the secular equation 

gyfoi) _ Ri rA8 , 

Again, the coefficients h 2 and /14 arc then found imposing 
the same boundary conditions for a fixed energy. 

Appendix B: Meaning of spin in topological 
insulators 

The bands of a semiconductor, in the presence of spin- 
orbit coupling, are characterized by their real total angu- 
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(a) (b) 




r (nm) r (nm) 

FIG. 10: (Color online) Radial probability associated with the first and third mode for a ring with a — 1000 nm and w = 50 nm 
and with negative Dirac rest mass M — — 10 meV. 



lar momentum (meaning orbital plus real spin). For our 
topological insulator we can express the relevant bands 
at the T-point as^ 



r *4 
r *4 



r 7 ,±- 



i 

±2 



T2 {X±lY) 



±1 
2 



± 



1 



V6 
1 



(X ± iY) 
(X±iY 



T~ 2 )T2Z 



1 

±- 

2 



1 



±- 



where each one of the four pairs obeys = 

— T|ri,+|), with T the time-reversal operator. The 
states considered in the effective 4-band model are de- 
scribed as^i 



We note that the two blocks r = ±1, although often re- 
ferred to as spin up and spin down blocks, do not exactly 
correspond to eigenstates up and down of the real spin 
along the Z axis. We can indeed, using Eq. (|B1|) . calcu- 
late the real spin expectation value of these two states, 
obtaining 



+ -A (-***! + *** 2 )f + 

* 2 | 2 [i + H 2 ] + ^)^ 

-\S\-) = -(+|!|+). (B4) 



\Ei±) 
\Hi±) 



a 



T 6l ±-)+f5 



r -4 



r 8 ,±- 



with \E\±) and \H\±) two set of Kramer's partners. We 
can therefore explicitly express the spinor on the block 
t = 1 as 



*a|Hi- 



and similarly we can express its time reversal (we use 
here the unsymmetrized basis) 



\-) = n+) = -n\Ei-) - n\Hi-) 



It is now clear that if we choose a spin quantization 

axis along the orthogonal |±) states will corre- 

(Bl) spond to spin up and spin down eigenstates, respectively. 

Therefore, in principle, particles injected with spin along 
this direction will enter only |+) states, conversely the 
ones with opposite spin will enter only |— ) states. How- 
ever the spin quantization direction depends implicitly 
on energy through "Fi and \&2- In conclusion, in a small 
(B2) enough energy range we can consider |±) as spin states 
along the spin quantization axis given by Eq. (|B4|) and 
it is meaningful to think of spin injection or detection 
processes, where the TI ring is coupled to external spin- 
selective materials like ferromagnetic contacts or helical 
(B3) edge states of a two dimensional TI. 
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